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THE  DOWLING  WALL  PRESSURE-SPECTRUM  ANALOGY  APPLIED  TO  AN 
ISOTROPIC  TWO- LAYERED  ELASTIC  MEDIUM 

Introduction 


The  relative  turbulent-boundary- layer  wall  pressure-spectrum  character¬ 
istics  on  an  elastic  medium  were  predicted  by  Dowling  through  the  solution  of 
the  Lighthill  equation  using  Green's  functions.  The  analogy  is  extended  in 
this  report  to  include  a  general  two-layered  isotropic  elastic  wall.  The  elas¬ 
tic  wall  is  treated  as  a  viscoelastic  material  through  the  use  of  complex 
shear  and  Young's  moduli.  The  relative  pressure-spectrum  characteristics  are 
given  on  the  top  surface  of  the  layer  of  the  elastic  medium  in  contact  with 
water  and  also  at  the  surface  between  the  two  elastic  layers.  Air  is  assumed 
to  be  the  fluid  in  contact  with  the  back  surface  of  the  two- layered  elastic 
wall.  The  effect  of  elastic-wall  moduli  and  damping  on  the  strength  and  loca¬ 
tion  of  surface-pressure  singularities  is  shown  to  be  rather  significant. 


BACKGROUND 


.0  r  4.1 


i- 


The  generation  of  noise  produced  by  turbulent  flows  over  elastic  surfaces 
is  the  subject  of  much  current  research.  Flow  noise  over  aircraft  bodies,  jet 
noise,  and  noise  generated  by  flow  over  underwater  vehicles  are  but  a  few 
examples  of  problems  in  this  area.  In  this  report,  we  have  extended  the  work 
of  Dowlingl-'5  to  the  analysis  of  the  relative  wall  pressure  spectrum  on  the 
surface  of  a  two- layered  elastic  medium  under  a  turbulent  boundary  layer  in 
water.  The  other  surface  of  the  elastic  medium  is  exposed  to  air  at  rest. 

Such  an  infinite  configuration  approximates  that  which  occurs  in  a  water  tun¬ 
nel  with  isotropic  structural  walls  covered  by  an  isotropic  viscoelastic  mate¬ 
rial.  In  the  next  section  of  this  report,  we  derive  the  wall  pressure  spectrum 
due  to  turbulent  flow  by  means  of  a  Green-function  analogy  applied  to  Light- 
hill's  equation  describing  the  density  field  in  the  fluid.  The  following  sec¬ 
tion  incorporates  the  properties  of  the  elastic  layers  into  the  expression  for 
the  surface  pressure  spectrum  through  the  boundary  conditions  at  the  interface 
between  the  fluid  and  elastic  medium.  Our  model  also  includes  the  effect  of 
turbulent-boundary- layer  shear  stresses  at  the  elastic  surface.  In  the  final 
section,  we  examine  the  structure  of  the  pressure  spectrum  on  a  bare  steel 
plate  at  5  and  10  kHz.  Inen,  we  compare  this  to  the  pressure  spectrum  on  the 
surface  of  two  different  viscoelastic  layers  applied  to  the  steel  plates.  The 
results  show  that  the  propagation  speeds  of  compressive,  shear,  and  flexural 
waves  in  the  elastic  layers  govern  the  location  of  singularities  in  the  pres¬ 
sure  spectra.  It  is  also  shown  how  damping  can  affect  the  strength  of  the 
singularities  and  the  level  of  the  local  maxima.  Q  ^ _  \  ^ 


-”■*> 
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SOLUTION  TO  LIGHTHILL  EQUATION 


The  analogy  derived  for  predicting  the  pressure  spectrum  on  a  single 
elastic  layer  through  the  solution  to  the  Lighthill  equation  and  elasticity 
equations  is  given  in  detail  by  Dowling. 3  In  this  section,  we  outline  the  gen 
eral  method  of  solving  the  Lighthill  equation  so  we  can  show  how  the  pressure 
spectrum  couples  to  the  two-layered  elastic  medium  in  the  following  section. 


Our  model  is  shown  in  figure  1.*  For  the  flow  speeds  of  interest  in  our 
underwater  applications,  the  effect  of  the  convected  flow  field  is  small  and, 
therefore,  does  not  appear  as  a  parameter  in  the  final  analogy.  Likewise,  the 
boundary- layer  thickness  appears  as  a  nondimens ional  parameter  in  the  final 
equation  for  the  pressure  spectrum  and  not  as  an  explicit  variable.  Also,  we 
consider  only  surfaces  whose  mean  displacement  is  linearly  disturbed  from  rest 
position  as  this  allows  us  to  make  some  simplifications  in  the  application  of 
the  Green- function  solution  at  the  surface  and  in  linearizing  the  momentum 
equation.  The  assumptions  and  simplifications  above  allow  us  to  proceed  to 
solve  the  relative  turbulent-boundary-layer  wall-pressure  problem. 


The  Lighthill  equation  describing  the  model  in  figure  1  is 

32 


C4  -  «*»*) 


-  C2V2)(p  -  p„) 


ayl3yj-{pvi,i 

♦  [p  - 


(1) 


Considering  the  source  term  on  the  right-hand  side  as  a  linear  disturbance,  we 
can  apply  the  technique  of  Green  functions  to  solve 


-  c2V2^G(y ,t  jx,t)  »  6(2  -  ?,t  -  t)  .  (2) 

Superposition  of  this  Green  function  with  the  source  term  will  yield  the  den¬ 
sity  perturbations  in  the  fluid  and  on  the  surface  of  the  upper  elastic  layer. 


The  density  perturbation  given  by  taking  the  integral  over  the  entire 
volume  yields 


(p  -  p„)(x,t) 


ay^yj 


d3ydt  . 


(3) 


Application  of  Green's  theorem  to  equation  (3)  will  give  the  volume  and  surface 
integrals. 


•All  figures  have  been  placed  together  at  the  end  of  this  report. 
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(P  -  Pg) 


32G 

3yi3yj 


3G 

0V2  3t 


d3ydt 


a 


3G 

2a  3ya 


♦ 


(4) 


The  quantities  of  particular  interest  to  us  are  the  Fourier  transforms  of 
the  pressure  spectrum  taken  in  time  and  in  the  two  horizontal  directions.  The 
Fourier  transform  taken  on  the  surface  integral ,  only  (the  overbars  denote 
transformed  variables),  is  shown  by 


(p 


P0)(x.t) 


t 


ij 


32G 

*1*1 


dydr  ♦ 


p(0,k,u>)^-(0,-k,-w) 


+  i£p0uv2(0,k,w)  ♦  ka52a(0,k,(j)j  6(0,-k,-u))|d2kdu)  .  (5) 

If  the  Green  function  is  chosen  to  make  the  surface  integral  vanish,  we 
arrive  at  a  boundary  condition  for  this  Green  function,  which  is 

^■(y2,-k,-u»)  ♦  iN (k , w) <5 (y2 , -k , -w)  ■  0.  at  y2  «  0  ,  (6) 

where  N (iT.w)  *  p  m2  +  k  Y  .  The  functions  Z  and  Y  represent  the  generalized 

^  ^ 

impedance  functions,  Z(k,w)  *  v2(0,k,u)/p(0,k,u))  and  Ya  =  a2a(0,k,u))/p(0,k,w) . 
These  impedance  functions  imply  that  the  elastic  deformation  and  fluid  motion 
can  be  described  in  this  linear  fashion. 

The  Green  function  satisfying  the  boundary  condition  in  equation  (6)  can 
be  obtained  directly  and  its  derivative,  32/3y^3yj,  is  given  as 


32G 

3yi3yj 


(y ,t | x,t) 


-1  -  n_ii  e-i«-(t-t)-ika(ya-xa)-iYy2d2i.dw 

(2ff)3ic 2J  Y-N 


(7) 


The  function  D^j  is  given  as 


°l,2 

°2,1 

ft 

tt 

°1.1  ■  K! 

°2,2 

°3,1 

"  °1,3 

*  M3 

°3,2  “  °2,3  *  K3Y 

°3,3 

(8) 


The  function  y  is  defined  as  ^«2/c2  -  K2  -  K2)1/2  and  the  appropriate  roots  to 

choose  for  y  when  it  is  real  or  imaginary  depend  on  the  choice  of  branch  cuts 
used  in  evaluating  the  Green  function. 
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Inclusion  of  equation  (7)  into  equation  (S)  will  yield,  after  taking  the 
transform, 


P(k,u)  =  - 


i(y  -  N) 


/■ 


e'1Yy2Tij(y2,k,w)dy2 


(9) 


In  equation  (9) ,  the  pressure  perturbation  has  been  related  to  the  density 
change  through  the  expression  p  *  c2(p  -  pQ).  This  relation  implies  a  simple 

linear  compressible  (i.e.,  acoustic)  relation  of  pressure  fluctuations  to  den¬ 
sity  fluctuations. 

To  solve  for  the  total-pressure  spectrum,  we  take  the  cross  correlation 
of  both  sides  and  obtain 


P(k,u.) 


Dij°U 


lYy2+1Y  y2Tijjtt(y2 ,y2 ,k,ui)dy2dy2 


(10) 


The  nature  of  Tjj^,  the  Fourier  transform  of  the  turbulent  source  terms, 
is  believed  to  be  such  that,  in  the  low-wavenumber  region  of  interest,  the 
effect  of  the  small  displacements  at  the  elastic  surface  and  fluid  compressi¬ 
bility  do  not  appreciably  alter  the  spectrum  of  the  turbulent  source  terms. 
Therefore,  one  can  nondimens ionalize  the  pressure  spectrum  to  obtain  the  final 
expression  as 


P(iU) 


PijPU 
|y  -  N | 2 


p2U3h5Q 


ijkZ  * 


(ID 


where  h  is  the  boundary-layer  thickness  and  U  is  the  free  stream  velocity 
above  the  boundary  layer.  The  superscript  plus  denotes  the  complex  conjugate 
and  is  the  Fourier  transform  of  the  cross  correlation  of  the  turbulent 

source  terms.  To  obtain  the  nature  of  the  wall  pressure  spectrum,  it  remains 
only  to  relate  the  elastic  properties  to  the  impedance  functions  Z  and  Y  . 


ELASTIC  IMPEDANCE  FUNCTIONS 


The  elastic  media  described  in  this  model  are  homogeneous,  isotropic,  and 
infinite  in  the  horizontal  directions.  Therefore,  we  can  simplify  our  model 
by  considering  one  horizontal  component  only,  since  the  other  horizontal  com¬ 
ponent  can  be  obtained  by  a  coordinate  transformation  applied  to  the  reduced 
problem.  We  are  essentially  working  in  two-dimensional  space. 

For  elastic  waves  in  layered  media,  the  velocities  are  given  by 
Brekhovskikh.4  Since  v  *  iwu,  the  displacements  will  be  given  here  as 
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u  *  7$  7xi£  , 

where  <t>  and  are  functions  of  the  material  properties  of  the  media; 


=  A  si 


"($ '  k')'/2y  * 8  C0S(?F '  k>) 


1/2 


and 


=  C  sin 


‘"(f '  kf)'/2y  * "  cosfe ' k?) 


1/2 


(12) 


where  Cj  and  are  the  compressional  and  shear  wave  speeds  in  the  i-th  layer. 
A,  B,  C,  and  D  are  as  yet  undetermined  coefficients  that  will  be  resolved 
through  application  of  the  boundary  conditions.  The  wave  speeds  can  be  related 
to  the  elastic  material  constants. 


where  and  are  the  Lame  constants  corresponding  to  the  i-th  layers  and 
is  the  mass  density  corresponding  to  the  i-th  layer. 

With  the  two-dimensional  situation,  there  are  a  total  of  four  boundary 
conditions  to  be  applied  to  each  interface,  continuity  of  normal  stress,  tan¬ 
gential  stress,  normal  displacement,  and  tangential  displacement.  Denoting 
the  upper  fluid  properties  with  no  superscript  and  the  upper  medium,  lower 
medium,  and  lower  fluid  with  superscript  prime,  double  prime,  and  triple  prime, 
respectively,  we  arrive  at  12  boundary  conditions. 


BOUNDARY  CONDITIONS 


At  the  Surface  y 


Boundary  conditions  at  the  surface  y  =  0  are 

*2 


■F  =  o'  =  A’fi-1-  +  - 

22  \  3x2  3y2 


-  [A*(kf  *  y’2)  ^  2u ' y ' 2] B '  ♦  2iu'k^C'  , 
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a  =  a ' 
12  12 


U'2(^  +  Tt  '  ifi)  "  2ilJ,kly,A,  +  W'^2  *  kt)D' 

V2  =  V2  =  iw(l7  +  iSf)  ’  iuY'A’  -  u^D’  , 


ri s  vi 3  iw(^f  ■  i't)  =  -wkiB’  •  iuYsc'  • 


At  the  Surface  y  a  T \ 

Boundary  conditions  at  the  surface  y  *  Tj  are 

°22  *  °22*  ”[*'(ki  *  Y'2)  +  2u,Y,2j(A'  sin  y'T1  ♦  B'  cos  y'T^ 


2iw'k1y;| 

lc' 

cos  y'T.  -  D* 
s  1 

sin 

-  [v,(k; 

♦ 

y"2^  ♦  2u"Y"2j 

(a"  : 

2iu"k  y"i 
1  s 

(c" 

COS  y"T  -  D" 
s  l 

sin 

a12  *  < 

J12 

;  2iu 'k^y' (a* 

cos 

♦  1 

Yg2  -  k2)  (c  sin  y 

=  2iu"k1Y"(A"  cos  y"Tj  -  B"  sin  y'Tj 
♦  u"(y's'2  -  k2)  (C"  sin  Y'g'Tj  +  D"  cos  yJJtJ  , 

*  v” ;  y'(a*  cos  y'T1  -  B'  sin  y'T^  ♦  ik^C'  sin  Vg^  +  D'  cos  YgT 
=  y"(a"  cos  y,rTl  -  B"  sin  y"T^  +  iki(c”  sin  y's'T1  +  D"  cos  y'^Tj  , 
i 

v |  =  v”;  ik^A'  sin  y'T1  ♦  B'  cos  y'T^  -  Yg(c  cos  y^  ♦  D'  sin  y^T 


=  ik, (a" 


sin  y"T,  B"  cos  y"Tf)  -  y"(c"  cos  y’T. 


-  D"  sin  y"T  ) 
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At  the  Surface  y  *  T2 

Boundary  conditions  at  the  surface  y  =  T2  are 

-p'"  =  0»2  =  .  +  Y"2)  ♦  2w"y"2]  (am  sin  y"T2  +  B"  cos  y"T2) 

+  2iu"k,Y"(cM  cos  y”T,  -  D”  sin  y"T,^  > 

1  s  \  S  2  S  2/ 

a12  “  al2  =  2iw"k1Y"(A"  cos  y"T2  -  B"  sin  y"T2) 

+  u"(y's'2  -  ki)  (c"  r,'*n  y^'T2  +  D"  cos  y'^T2)  , 
v'2"  =  v^’  =  iu)Y"(A"  cos  y”T2  -  &"  sin  y"T2)  -  wk^C"  sin  y^  +  D"  cos  y”T_,)  , 


v'"  =  v"  =  -uk  f A"  sin  y"T  -  B"  cos  y"T„\  -  wfc"  cos  y"T,  -  D"  sin  y'^O  . 
11  l\  2  2/  S\  S  2  's  2/ 


FURTHER  REDUCTION  OF  THE  SYSTEM 

Linearization  of  the  momentum  equation  in  the  fluid  results  in  two  rela¬ 
tionships  that  further  reduce  the  system  of  equations.  These  relationships 
express  tangential  velocity,  normal  stress,  and  tangential  stress  as 

-ol2  kjp 

v,  *  - , - .  (14 

1  p(iujv)l/2  0W 

On  the  lower  surface,  since  the  lower  fluid  is  air  at  rest,  the  relationship 
between  normal  velocity  and  normal  stress  can  be  described  as 

V,"  =  X_E_  .  (is 

2  p’"w 

Using  equations  (14)  and  (15)  in  the  boundary-condition  equations  and  assem¬ 
bling  the  result  in  matrix  form,  we  arrive  at 
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The  system  of  equations  can  be  solved  by  inverting  the  A  matrix  and  dividing 
the  unknown  variables  by  P,  the  surface  pressure  on  the  upper  elastic  layer, 


The  coefficients  of  A  are  given  as 

al,2  s  '[A'(kl  +  Y'2)  +  2u'y’2]  , 

ai,3  =  2i**v;  * 

a„  ,  =  2iu'k1Y'  , 
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■wk  j  » 

-1  , 

*u>k  j  f 

-iwYg  » 

l/p  (i<i."0  l/2  , 


-Hk?  * 

y*1) 

|  ♦  2u'y'2 

'jsin  y ,Tl  , 

-{>'(k5  * 

y*  0 

I  +  2u'y'^ 

■| cos  y'Tj  , 

2iii'kjYg 

COS  1 

'iTl  * 

sin 

^;ti  > 

[‘"H  * 

y..2) 

+  2u"y"2 

sin  y'Tj  , 

[»"(k5  * 

y"2) 

« 

+  2|i’’Y"2 

jcos  y"T1  , 

-2iu”klY's' 

cos 

ysTi  ' 

2i w”k iy*s'  sin  y^  , 
2iu'lc  y'  cos  y'T  , 
-2ip'k  y'  sin  y'T^  , 
y'(Y;  -  k*)sin  Y;Tl  . 

y,(Ys  ■  ki)C0S  >sTl  * 
-2in"kjY"  cos  y"Tj  , 

2iw"klY"  sin  y"Tj  , 

-u"(y;  -  kf)sin  y^T j  . 

-m"(yL'  -  k*)cos  y'T.  , 
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*10,6  *  '2iu,V"  Sin  Y"T2  ' 

*10,7  *  »"(y?  '  kl)sin  YsT2  * 

*10, a  u"(Y’s2  *  k?)cos  ysT2  • 

*10,11  m  -1  ' 

*11,5  *  iwY"  C0S  Y’*T2  > 
an,6  *  -iwy"  sin  y"T2  , 

*n,7  *  -“ki  sin  ^2  • 

*n, a  *  'uki  cos  ys'T2  • 

*11  12  *  -Y'"  p"'/pmW  , 

*12 , 5  *  *wkl  Sin  ^2  * 

*12,6  *  “kl  C0S  Y"T2  * 

*12,7  *  ” iwY  s  C0S  Ys'T2  ’ 

*12,8  *  iwYS  Sin  ysT2  • 

*12  11  *  l/0'"(iu>v'")l/2  , 

and 


*12,12 


k1/p,,,to  . 


All  unlisted  coefficients  are  zero.  Since  the  impedance  functions  are 
desired  to  complete  the  pressure  analogy,  we  obtain  them  from  the  matrix  equa¬ 
tions  and  express  them  in  terms  of  the  inverted  matrix  A, 


and  * 


t"  v>  ■  HO  H  •  ■ 

where  atj  indicates  that  particular  entry  in  the  inverse  matrix  A'1 
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Equation  (11)  for  the  pressure  spectrum  can  now  be  evaluated  by  substi¬ 
tuting  the  values  of  Z(k,u>)  and  Y(k,w),  as 


P(k,u>)  = 


DijDk* 


y  *  pu,(ai5,i) +  ki(aio,4)  *  ki(a;!i) ♦  ^ir(a;S)l 


02U3h5Qs 


ijkt 


(18) 


PRESSURE  SPECTRUM  RESULTS 


Since  the  manual  inversion  of  a  12  by  12  complex  matrix  is  a  formidable 
task,  the  solution  to  equation  (17)  was  performed  by  a  double-precision  com¬ 
plex  matrix-inversion  subroutine  from  the  commercially-available  international 
mathematical  and  statistical  libraries  (IMSL)  on  a  VAX  11/780  computer.  To 
check  the  accuracy  of  the  solution  algorithm,  the  two  cases  of  a  single  steel 
plate  of  1  or  5  cm  backed  by  water  at  10  kHz,  given  by  Dowling,  are  compared 
with  similar  cases  here.  The  pressure  spectrum  Dowling  shows  for  those  cases 
are  given  on  the  lower  face  of  the  elastic  slab.  To  obtain  a  similar  pressure 
spectrum,  we  must  alter  equation  (18). 

From  equation  (17),  the  ratio  p'”/p  is  given  by 


p'» 


(19) 


Substituting  equation  (19)  into  equation  (18)  gives  the  pressure  spectrum  on 
the  lower  surface, 


P(k,w) 


D. 

i: 

iDia  '  ( a i 2 , i ) 

-  —fa" 1  1 

pw\  12,4> 

CM 

Y  ♦  ow(a-3> 

i)  *  k i(aI o , *♦) 

*  k l(a9  >  l) 

p2U3h5Q 


ijkn 


(20) 


Figures  2  and  3  are  plots  of  a  nondimensional  function,  F,  which  shows 
the  characteristics  of  the  pressure  spectrum  on  the  lower  surface.  The  func¬ 
tion,  F,  is  defined  as 


PW  (ai2 , 4  )| 

Y  +  PW(alJ,l)  *  kl(al0,4) 

*  ki(s!i) 

*=W.)I 

(21) 


The  plots  show  20  log  F  versus  nondimensional  K-*-  *  kc/w.  These  plots  compare 
exactly  with  those  in  Dowling^  of  1-  and  5-cm  bare  steel  plates  backed  by  water 
at  5  and  10  kHz,  respectively. 

The  major  application  of  the  wall-pressure  analogy  in  this  report  will  be 
to  investigate  the  surface-pressure  spectrum  given  by  equation  (18)  for  cases 
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where  the  upper  fluid  is  water,  the  lower  fluid  is  air,  and  the  medium  is  a 
combination  of  viscoelastic  and  steel  materials.  The  pressure  spectrum  is 
given  both  on  the  top  surface  of  the  upper  layer  and  on  the  interface  between 
the  two  elastic  layers.  As  a  base-case  comparison  for  the  material  combina¬ 
tions  to  follow,  we  show  the  surface-pressure  spectrum  on  the  upper  surface  of 
a  1-in.  steel  plate  at  5  and  10  kHz  in  figures  4  and  5.  All  pressure  spectra 
plotted  in  this  report  are  similar  to  those  in  figures  2  and  3  in  that  they  are 
plots  of  the  nondimens ional  function  F  evaluated  at  the  appropriate  surface. 

To  plot  the  absolute  pressure  magnitudes,  we  would  need  to  know  the  values  of 
the  function.  Since  very  few  measurements,  if  any,  have  been  made  of  the 

turbulent-source  cross-correlation  terms,  our  intention  here  is  to  show  rela¬ 
tive  reductions  or  increases  in  pressure  spectrum  and  not  absolute  magnitude. 
Also,  the  structure  of  the  low-wavenumber  content  of  the  pressure  spectrum  can 
be  examined  for  local  maxima  and  minima. 

It  is  clear  from  figures  2  through  5  that  the  pressure  spectrum  always 
contains  a  discontinuity  at  K*  »  1.  At  5  and  10  kHz  in  figures  4  and  5, 
another  discontinuity  appears  moving  from  right  to  left  across  the  wavenumber 
spectrum  as  the  frequency  increases.  For  5  kHz,  the  flexural  wavelength  of  a 
1-in.  steel  plate  in  vacuo  corresponds  to  K+  *  1.36..  At  10  kHz,  this  value  is 
close  to  K*  a  1.  One  can  see  that  the  flexural  wavelength  of  the  plate  intro¬ 
duces  another  singularity  into  the  pressure  spectrum.  This  singularity  moves 
to  lower  wavenumbers  with  increasing  frequency  until  it  reaches  K+  *  1.  At 
higher  frequencies,  the  singularity  does  not  appear  to  move  below  K+  *  1.  Elas 
tic  surfaces  that  have  a  flexural  propagation  speed  greater  than  that  of  an 
acoustic  wave  in  the  fluid  do  not  appear  to  have  a  double  singularity  in  the 
pressure  spectrum  due  to  their  flexural  waves.  This  detailed  behavior  is  con¬ 
sistent  with  general  theoretical  observations  of  Dowling.2 

Now,  we  can  proceed  to  investigate  the  relative  dB  gains  or  reductions  in 
the  turbulent  wall  pressure  spectrum  as  a  result  of  viscoelastic  surfaces 
applied  to  the  1-in.  thick  bare  plates.  The  viscoelastic  properties  we  refer 
to  are  obtained  by  using  a  complex  shear  and  Young’s  modulus  with  a  specified 
thickness  and  mass  density.  Figures  6  and  7  show  the  top-  and  mid-surface 
pressure  spectrum  (mid  surface  is  the  interface  between  the  two  elastic  media) 
for  a  1-in.  steel  plate  at  10  kHz  covered  with  3  in.  of  viscoelastic  material 
having  the  properties  G  =  43  ♦  4.3i  psi,  E  *  129.1  ♦  12. 9i  psi,  and  p  *  0.0001 
lb-s2/in.4.  This  material  represents  a  fairly  "compliant"  surface  with  a  den¬ 
sity  near  that  of  water. 

The  singularity  due  to  the  flexural -plate  speed  is  now  seen  to  be  reduced 
in  magnitude  in  the  top-surface  spectrum  at  10  kHz  from  that  of  the  bare  steel 
plate  but  remains  quite  distinct  in  the  mid-surface  spectrum.  In  the  top- 
surface  spectrum,  the  elastomer-coated  plate  demonstrates  a  noticeable  reduc¬ 
tion  in  the  pressure  spectrum  for  wavenumbers  0.5  <  K*  <  1.  There  is  a 
slightly  larger  reduction  for  K*  >  1.1.  This  suggests  that  this  candidate 
material  may  actually  reduce  the  level  of  turbulent  "noise."  However,  in  fig¬ 
ure  7,  which  shows  the  mid-surface  pressure  spectrum,  the  dramatic  reduction 
from  that  on  the  bare  plate  corresponds,  for  the  most  part,  to  transmission 
loss  through  the  layer.  Also  appearing  in  figure  7  is  a  singularity  occurring 
at  K*  a  0.3.  This  singularity  actually  occurs  in  most  of  the  plots  and  corre¬ 
sponds  to  the  compress ive-wave  speed  of  the  steel  plate  given  by 
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Except  for  this  singularity,  the  supersonic  regime  of  the  top-surface 
pressure  in  the  figure  6  spectrum  appears  to  have  a  very  flat  profile.  The 
rise  in  pressure  level  with  movement  from  0  <  K*  <  1  for  the  bare  steel  plate 
does  not  occur  in  the  case  of  this  covered  plate.  This  particular  case  indi¬ 
cates  a  reduction  in  the  top-surface  pressure-spectrum  level  for  some  super¬ 
sonic  wavenumber  elements  and  for  all  the  subsonic  wavenumbers  shown  in  this 


Figures  8  and  9  show  the  same  3-in.  viscoelastic  coating  at  5  kHz.  Again, 
we  see  a  reduction  in  the  surface  pressure  at  the  top  surface  of  the  coating 
from  that  on  the  bare  steel  plate  and  a  significant  transmission  loss  is  evi¬ 
dent  in  the  mid-surface  pressure  spectrum.  The  elastomer  appears  to  have  less 
of  an  effect  on  the  pressure  spectrum  at  lower  frequencies  in  that  the  spec¬ 
trum  does  not  attenuate  as  dramatically  (on  the  mid-surface  at  higher  wave¬ 
number)  for  5  kHz  as  it  did  for  10  kHz.  On  the  top  surface,  the  flat  level  of 
the  pressure  spectrum  is,  again,  seen  to  be  in  the  supersonic  regime.  The 
examples  in  figures  6  through  9  show  the  frequency  dependence  of  the  pressure 
spectrum  for  a  fixed  set  of  material  properties.  It  is  more  likely,  though, 
that  real  material  properties  would  also  change  with  frequency. 

Figures  10  and  11  show  results  for  a  1-in.  steel  plate  coated  with  an 
elastic  material  layer  having  the  properties,  at  S  kHz,  G  *  90,000  psi,  E  = 
252,000  psi,  and  p  *  0.0001  lb-s2/in.4.  In  figures  12  and  13,  5  percent  damp¬ 
ing  has  been  added  to  the  above  shear  and  Young's  moduli.  The  purpose  of  using 
this  particular  material  is  to  demonstrate  the  effect  on  the  pressure  spectrum 
when  a  steel  plate  is  coated  with  a  material  having  its  flexural,  compressive, 
and  shear  speeds  all  in  the  vicinity  of  the  acoustic-wave  speed.  When  we  try 
to  find  surfaces  that  reduce  turbulent-noise  levels,  we  see  that  a  surface 
with  propagation  speeds  near  the  wavenumber  region  of  interest  will  introduce 
additional  singularities  in  the  pressure  spectrum,  as  is  demonstrated  by  this 
particular  material. 

One  of  the  most  noticeable  features  in  figures  10  through  13  is  the 
increase  in  the  pressure  spectrum  at  higher  wavenumbers  over  that  of  the  bare 
steel  plate.  This  occurs  for  the  damped,  as  well  as  the  undamped,  case.  In 
addition,  the  flat  level  of  the  pressure  spectrum  on  the  top  surface  for  super¬ 
sonic  wavenumbers,  seen  for  the  softer  viscoelastic  surface  shown  in  figures 
6  and  8,  is  no  longer  shown.  Instead,  the  pressure  reaches  a  local  maximum 
around  K*  *  0.75.  If  one  looks  at  the  undamped  top-surface  spectrum  versus 
the  undamped  mid-surface  spectrum,  figures  10  and  11,  it  is  evident  chat  the 
pressure  spectrum  is  greater  on  the  mid  surface  than  on  the  top  surface  for 
K*  <  0.75.  However,  with  movement  out  to  >  2,  the  pressure  spectra  at  both 
surfaces  are  nearly  equal  in  magnitude. 

The  effect  of  a  small  amount  of  damping  in  shear  and  Young's  moduli  is 
most  dramatic  in  figure  13,  the  mid-surface  pressure.  The  structure  of  the 
pressure  spectrum  for  1  <  H-1,  <  1.75  has  changed  dramatically  from  the  undamped 
case  in  figure  11.  Nulls  not  seen  before  in  the  pressure  spectrum  are  now 
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evident.  The  nulls  probably  also  exist  in  the  undamped  case.  However,  the 
very  narrow  extent  of  these  nulls  is  such  that,  for  the  undamped  case,  the 
finite  points  used  to  plot  the  curve  are  spaced  too  far  apart  to  reveal  the 
nulls.  Damping  serves  to  broaden  the  wavenumber  extent  of  the  nulls,  giving 
enough  span  so  that  the  nulls  now  appear. 

In  both  the  mid-  and  top-surface  pressure  spectra,  damping  in  the  coating 
appears  to  control  the  magnitudes  of  the  singularities  but  does  not  change 
their  location  on  the  wavenumber  axis.  It  would  appear  that  very  heavy  damp¬ 
ing  might  remove  the  singularities  for  K+  >  1.  The  surface-pressure  spectrum 
for  the  lightly  damped  coating  of  figure  12,  when  compared  to  that  for  the  bare 
steel  plate,  shows  a  marked  difference.  However,  it  is  difficult  to  assess 
visually  whether  the  coating  results  in  more  or  less  turbulent  noise  because 
of  the  complexity  of  the  wavenumber  spectra.  One  final  interesting  feature  is 
the  fact  that  damping  increases  the  level  of  the  local  maximum  in  the  pressure 
spectrum  at  K+  =  0.75  for  both  the  top  and  mid  surface.  It  is  not  certain  why 
this  is  so. 


CONCLUSIONS 


The  Dowling  pressure  analogy  for  an  elastic  medium  under  a  turbulent¬ 
boundary-layer  flow  has  been  extended  to  a  two-layered  elastic  medium.  This 
analogy  predicts  the  nature  of  the  relative  pressure  spectrum  at  any  surface 
in  the  elastic  medium  and  can  demonstrate  graphically  the  results  of  different 
combinations  of  isotropic  elastic  and  viscoelastic  materials.  The  assumption 
applied  in  the  analogy  to  the  turbulent-source  terms  limits  the  range  of  appli 
cability  of  the  model  to  low  wavenumbers  relative  to  the  convective  peak.  In 
addition,  large  surface  motion  would  violate  the  assumption  of  a  linear  rela¬ 
tionship  between  structural  displacement  and  fluid  velocity.  For  the  low- 
wavenumber  region,  this  model  provides  a  unique  way  of  looking  at  multilayered 
systems  and  their  behavior  relative  to  an  unspecified  turbulent-boundary- layer 
loading. 

Specifically,  we  have  shown  the  pressure  spectra  of  two  different  visco¬ 
elastic  materials  on  a  steel  plate.  The  first  viscoelastic  coating  chosen 
showed  local  reductions  in  the  pressure  spectrum  on  the  top  surface  and  on  the 
mid  surface.  The  second  viscoelastic  coating  was  chosen  to  demonstrate  that 
surfaces  which  have  propagation  speeds  near  the  acoustic-wave  speed  can  intro¬ 
duce  new  structure  into  the  pressure  wavenumber  spectrum,  similar  to  an  obser¬ 
vation  by  Dowling?  on  coatings  with  sound  speeds  less  than  that  of  the 
turbulent  fluid. 

It  was  also  shown  that  plate  flexural  and  compressional  singularities  can 
be  affected  by  damping  in  the  elastic  layers.  This  model  is  intended  to  pro¬ 
vide  a  means  by  which  the  performance  of  coatings  applied  to  parent  material 
can  be  investigated  as  to  their  influence  on  turbulent-boundary- layer  noise 
levels.  The  inclusion  of  shear  waves  in  the  elastic  layers  and  surface  shears 
due  to  the  turbulent  flow  make  the  model  more  complete  for  use  in  predicting 
wall  pressure  spectra. 
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Figure  1.  Geometry  of  the  Two- Layered  Medium  and  Flow  Field 
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